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Using a data sample collected with the BESIII detector operating at the BEPCII storage ring, 
we observe a new neutral state ^c(3900)'’ with a significance of 10.4cr. The mass and width are 
measured to be 3894.8 ± 2.3 ± 3.2 MeV/c^ and 29.6 ± 8.2 ± 8.2 MeV, respectively, where the first 
error is statistical and the second systematic. The Born cross section for e^e“ —>■ J/i/) and the 

fraction of it attributable to 7r°2c(3900)° —>■ in the range Ecm = 4.19 — 4.42 GeV are also 

determined. We interpret this state as the neutral partner of the four-quark candidate ^c(3900)*. 

PACS numbers: 14.40.Rt, 14.40.Pq, 13.66.Be 


A new charged charmonium-like particle Zc(3900)^ has recently been observed through its decay to 
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TT^J/tp by BESIII, Belle and a Northwestern Univer¬ 
sity group using CLEO-c data [H-Q. This state lies 
just above the threshold for DD production, simi¬ 
lar to the bottomonium-like resonances Zb(10610)^ and 
Zb(10650)^ that have been observed by Belle at an en¬ 
ergy just above BB* threshold 0. BESIII also ob¬ 
served a structure, Zc(3885)^, in the process e+e“ —>• 
, with mass close to Zc(3900)^@. Because 
the Zf- couples to charmonium and has electric charge, 
it can not be a conventional qq meson, but must contain 
at least two light quarks in addition to a cc pair. Pro¬ 
posed interpretations for include hadronic molecules, 
hadro-quarkonia, tetraquark states, and kinematic effects 
0 . The precise structures of the Zf- and other “XT Z" 
states remains unknown, and hence that their further 
study will lead to a deeper understanding of the strong 
interaction in the non-perturbative regime. 


Progress in clarifying this picture requires measure¬ 
ments of improved precision and searches for additional 
states. The first definitive observation of a neutral 
Zc state was a BESIII measurement of Zc(4020)° ^ 
7r°/ic i- Previously, 3 . 5(7 evidence for a candidate 
state Zc(3900)° decaying to Tr^J/ip was observed in re¬ 
port 0. In this Letter, we report the observation of 
Zc(3900)° in the process e+e“ —>■ 7 r° 7 r°J/^/> based on 
data collected with the BESIII detector at the BEPCII 
electron-positron collider. First measurements of the 
Born cross section for e“''e“ 7 r° 7 r° J/'0 and of the frac¬ 

tion of 7 r*’ 7 r°J /'0 production attributable to Zc(3900)° 
as a function of center-of-mass energy (Ecm) are also 
presented. Our data sample has an integrated lumi¬ 
nosity of 2809.4 pb“^ distributed over the Ecm range 
from 4.190 to 4.420 GeV Q, with an overall measure¬ 
ment uncertainty of 1.0%. The three largest samples have 
Ecm = 4.230 GeV (1091.7 pb-^), 4.260 GeV (825.7 pb-^) 
and 4.360 GeV (539.8 pb“^), with the remainder dis¬ 
tributed comparably among seven other energies Q. 

BESIII is a general-purpose magnetic spectrometer 
M with a helium-gas-based drift chamber (MDC), a 
plastic scintillator time-of-flight system (TOE), and a 
GsI(Tl) Electromagnetic Calorimeter (EMC) enclosed in 
a superconducting solenoidal magnet providing a 1.0 T 
held. The solenoid is supported by an octagonal hux- 
return yoke with resistive-plate counters interleaved with 
steel for muon identihcation (MUG). 


To study the signal response in the BESIII detector, we 
use a Monte Carlo (MG) package based on GEANT4 111 
to produce simulated samples for e“''e“ —>■ Z^ 

TT^J/ip and e+e" ^ without an intermediate 

zl in both cases with J/ip ^ e+e or . We gen¬ 
erate e+e” and Z^ n^J/ip with isotropic 

angular distributions. We simulate e+e” 7 r° 7 r°J /'0 
with a generator of hnal states with a J/ip and two 
pseudoscalars in EVTGEN IJ, [l3| and no intermedi¬ 
ate resonances contributing to the 7 r° production. 


To determine the Z^ mass resolution, a signal sample 
is generated at Ecm = 4.260 GeV with a Z^ mass of 
3.9 GeV/c^ and zero width. In measuring the cross sec¬ 
tion cr(e“''e“ 7 r° 7 r° J/'0) and ratio 

(7(e“*'e“ 7r°Zc(3900)° 7r°7r° J/'0) 

cr(e+e- ^ ttOttOJ/Pj) ’ 

MG samples for e+e” —>■ 7 r° 7 r°J/' 0 , with and without 
an intermediate Z^, and using mass and width values 
obtained in this analysis, are generated at all ten Ecm 
points. QED radiative corrections for J/ip —)■ £+£“ are 
incorporated with PHOTOS [1^, and initial-state radia¬ 
tion (ISR) is simulated with KKMC [l^ using the same 
parameters as in Ref. [i| . To study background, a generic 
Y (4260) sample and a sample of simulated continuum qq 
production at Ecm = 4.260 GeV equivalent to 500 pb“^ 
are used, as in Ref. [ij. 

Charged tracks are reconstructed from MDC hits. To 
optimize the momentum measurement, we restrict the 
angular range of tracks to be |cos0| < 0.93, where 

9 is the polar angle with respect to the positron beam. 
We require tracks to pass within ±10 cm of the inter¬ 
action point in the beam direction and within 1 cm in 
the plane perpendicular to the beam. Electromagnetic 
showers are reconstructed by clustering EMC energy de¬ 
posits. Efficiency and energy resolution are improved 
by including energy deposited in nearby TOE counters. 
Photons are selected by requiring showers with minimum 
energies of 25 MeV for |cos0| < 0.8 or 50 MeV for 
0.86 < |cos0| < 0.92. The angle between the shower 
direction and the extrapolation of any track to the EMC 
must be greater than 5°. A requirement on the EMC tim¬ 
ing suppresses electronic noise and deposits unrelated to 
the event. Candidates for 7 r° —>• 77 decays are selected by 
requiring the diphoton invariant mass to be in the range 
100 < < 160 MeV/c^. 

We search for e+e" 7 r° 7 r° J/'0 in events with exactly 
two good oppositely charged tracks and at least four good 
photons. In reconstructing J/ip ^ e+e", electron candi¬ 
dates must satisfy E/p > 0.7, where E is the EMC energy 
and p is the momentum measured in the MDC. To sup¬ 
press the small two-photon and Bhabha background, the 
two-track opening angle is required to be less than 175° 
for any (e“) with cos0 > 0.5 {cosO < —0.5). In se¬ 
lecting J/ip —> we require both muon candidates 

to satisfy E/p < 0.3 and at least one to have associated 
hits in more than six MUC layers. 

We reconstruct J/tp candidates if the dilepton in¬ 

variant mass is within the J/'0 signal region (2.95 < 
Mu < 3.2 GeV/c^). We loop over 7 r° candidates and 
select the two that do not share photons and have the 
smallest = Xic + xlc’ where Xic 
X^ values for the two one-constraint (1C) kinematic fits 
to the TT® mass, and xic i® t^® tbe 4C fit to 

the 7 r® 7 r® J/'0 hypothesis requiring 4-momentum conser¬ 
vation. To suppress combinatorial background we require 
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that there be fewer than two combinations meeting 
the tighter tt® criterion of 120 < < 150 MeV/c^. 

To search for and suppress non-Tr^Tr^ J/t/j events, the 
event is subjected to a 7C fit, adding mass constraints for 
both TT^s and the J/tp to 4-momentum conservation. To 
improve resolutions, for events with Xyc <230, the 7C- 
constrained momenta are used to construct and 

M^o^o. We verified that resonant structures in the 7r°7r° 
mass spectrum, such as /o(980), do not produce a peak 
in the distribution. Figure [T] shows the ir^J/'p in¬ 

variant mass distribution in data and the MC-determined 
background for Ecm = 4.260 GeV. Each event 

appears twice, once for each 7r°. Background processes 
are estimated by MC to contribute ^ 12% of selected 
events, dominated by XJ/p{X ^ and multi-pion 

final states. Because the location of the lower peak de¬ 
pends on Ecm while the higher peak remains fixed, we 
interpret the excess near 3.9 GeV/c^ as ^^(3900)° pro¬ 
duction and that near 3.4 GeV/c^ as its kinematic reflec¬ 
tion. 



(GeV/c?) 


FIG. 1. Invariant mass distribution for tP Jjp candidates in 
Ecm = 4.260 GeV data (points). The dashed histogram shows 
the MG background and the solid histogram is the sum of this 
background and iPiP Jjp production not from Z'^. 


We extract the yields and resonance parameters of 
Ec(3900)*^ by performing an unbinned maximum likeli¬ 
hood fit simultaneously to the n^J/p mass distributions 
for the three high-statistics samples. The fit lower limit 
is set to 3.65 GeV/c^ to avoid double-counting. The sig¬ 
nal shape is an S-wave Breit-Wigner with phase-space 
factor pq, where p is the momentum in the 
frame and q is the J/p momentum in the Z^ frame. It is 
convolved with a resolution function consisting of three 
Gaussians with parameters set by fitting the zero-width 
e+e“ —>■ 7r°Zc MG sample at Ecm = 4.260 GeV (aver¬ 
age resolution « 6 MeV/c^). The background shape is 
an ARGUS function [l^. We use the same Breit-Wigner 
and resolution functions for all energy points because res¬ 
olution dependence on Ecm is determined by MG simula¬ 
tion to be very small. The ARGUS parameters are varied 


independently in the fit, except that the cut-off is based 
on Ecm- 

Figure [2] shows the simultaneous fit to the three n^J/p 
invariant mass distributions, which returns a Ec(3900)*^ 
signal with a statistical significance of 10.4(7 and a of 
176 for 151 degrees of freedom. Yields at Ecm = 4.230, 
4.260 and 4.360 GeV are 225.3±41.0, 83.2±20.5, and 
47.5±12.7, respectively, with a sum of 356.0±47.6. The 
Ec(3900)° mass and width values with statistical errors 
are 3894.8±2.3 MeV/c^ and 29.6±8.2 MeV, respectively. 



FIG. 2. The simultaneously fitted iP Jjp mass spectra (55 
bins in M^oj/P) for (a) Ecm = 4.230 GeV, (b) Ecm = 
4.260 GeV, and (c) Ecm = 4.360 GeV. Dots represent the 
data, solid lines represent the fitted results and dashed lines 
represent fitted backgrounds. 

We determine the cross section ratio R and the 
e+e” —?► TT^TT^J/p Born cross section as functions of 
Ecm by measuring yields of Z° {N{Z^)) and 7r°7r°J/'0 
{NTT^ J/p)). N{Z^) is determined with a simulta¬ 
neous fit of the TT^J/p mass spectra for all ten Ecm 
samples. The signal function is the same as for the fit 
to the high-statistics samples, with the Ec(3900)° mass 
and width fixed to the results of that fit. Background 
shapes are ARGUS functions with the cut-off based on 
Ecm and other parameters constrained to be the same for 
all points. 

To obtain NJ/p), the dilepton mass spectra for 
all energies are fitted simultaneously. The small peak¬ 
ing background from XJ/p{X p is treated as a 

systematic error. For this determination the 7C kine¬ 
matic fit including J/p mass constraints is inappropri- 
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ate and the 4C fit results are used. Events are selected 
with a cut of xic ^ ^0 based on an optimization con¬ 
sidering statistical and systematic uncertainties. Each 
signal shape is a Breit-Wigner convolved with a double- 
Gaussian. The Breit-Wigner is fixed to the width of the 
J/'f/’ and the mass is allowed to vary to allow for possible 
mis-calibration of the momentum scale for reconstructed 
tracks. The mean of the first Gaussian of the resolution 
function is fixed to zero, while the other parameters are 
varied. The background shape is a first-order Cheby- 
shev polynomial with free parameters. In this fit, the 
parameters of the double-Gaussian and the polynomial 
are constrained to be same for all energy points, except 
for the normalization factor. 

The fraction of 7r°7r°J/'0 production attributable to 
Zc(3900)° is determined with Eq. O where e(Z^) is the 
efficiency for extracting the signal by the fit to the 
TT^ J/'0 invariant mass distribution, and J/iji) and 

e2(7r°7r°J/'0) are efficiencies for determining 
yields by fits to dilepton mass distributions for processes 
without and with an intermediate respectively. 


R = 


NiZ°) 

<Z°c) 

<Z0) 



+ (iV(7r°7r° J/'i/i) — 


£2 (7r° 7r° J/V')) / Cl (7r° 7r° J/V') 


( 2 ) 




The observed cross section for e+e“ —>■ Tr^ir^J/ip is cal¬ 
culated using Eq. (U where C is the integrated luminosity 
and e(7r°7r° is the weighted average of the efficiencies 
for events with a Z^ (e2(7r°7r°J/'i/;)) and without a Z^ 
(ei(7r°7r°J/'0)). The branching ratios —>■ e+e“) 

and B{J/ip ^~) are taken from the PDG [I 3 . 


(Jobs = iV( 7 r° 7 r°J/^) j 


C X e{TT°TT°J/'ljj) X 


(6(J/V' ^ e+e”) -k B{J/ip 


( 3 ) 


The Born cross section is calculated with o'Born = 
o'obs/[(l -|- ^)(1 -I- (5™'^)], where (1 + ^) is a radiative cor¬ 
rection factor obtained with KKMC [l^ and (1 -f <5™'^) 
is a vacuum polarization factor following Ref. [l^. Note 
that due to initial state radiation to e^e“ resonant struc¬ 
tures such as F(4260), (1 -I- S) depends on Ecm- The 
inputs and results are listed in Table HI In cases where 
there is no statistically significant signal, the upper lim¬ 
its at 90% confidence level are provided. For N(Z^) and 
the errors and upper limits are statistical 
only. A cap of 1 is set on R values. Figure |3lja) and (b) 
show R and O'Born as functions of Ecm with error bars 
that are statistical only. 

We consider several sources of systematic uncertainty 
in the Ec(3900)° mass and width measurements. For 
the mass determination, the largest uncertainty is that 


FIG. 3. (a) R (see text) and (b) aBorn{e'^e Tv^TT^J/ip) as 
functions of Ecm. Error bars are statistical only. 


associated with the absolute track momentum scale, esti¬ 
mated to be 2.0 MeV/c^ based on the difference between 
the dilepton mass determined by the fit and the nominal 
J/'0 mass. Uncertainty due to the knowledge of the beam 
energy is estimated to be 1.7 MeV/c^ based on a study of 
e+e“ ^ ■ Adjusting the cut on Xrc ±30 changes 

the mass by 1.2 MeV/c^, which we assign as the system¬ 
atic uncertainty associated with the kinematic fit. To as¬ 
sess the uncertainty from the signal parameterization we 
change the phase-space factor from pq to (S-wave to 
P-wave) and find a 1.1 MeV/c^ change in the mass. Ad¬ 
ditional systematic effects associated with fitting-range 
dependence (0.8 MeV/c^), background-shape sensitivity 
(0.3 MeV/c^) and Ecm dependence (0.2 MeV/c^) con¬ 
tribute at a lower level, leading to an overall system¬ 
atic error in M(Ec(3900)°) of 3.2 MeV/c^. The mea¬ 
surement of r(Ec(3900)°) has a total systematic error 
of 8.2 MeV, which includes similarly sized contributions 
from the kinematic fitting procedure (4.6 MeV), back¬ 
ground shape (4.1 MeV), fitting range (3.9 MeV), and 
Ecm (3.3 MeV), with a smaller effect due to the signal 
model (1.7 MeV) and none from the absolute mass scale. 
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TABLE I. Efficiencies, yields, R = and ■K^'K^J/ij) Born cross sections at each energy point. For 

N{Z^) and N{-k^'k^ J/ip) errors and upper limits are statistical only. For R and cTBom, the first errors are statistical and second 
errors are systematic. The statistical uncertainties on the efficiencies are negligible. Upper limits of R (90% confidence level) 
include systematic errors. 


Ecm 

(GeV) 

£ 

( pb-b 

(%) 

€l{7r'^7T^J/lp) 

(%) 

e2{7T^'rT^J/'ip) 

(%) 

€{7Z^ J / Ip ) 

(%) 

w(zD 
( 90 % C . L .) 

C 

c 

R 

( 90 % C . L .) 

1±<5 

1 ± 5 ”“ 

tTBorn ( pb ) 

4.190 

43.1 

20.8 

20.4 

20.1 

20.2 

< 11.1 

8.2 ± 3.0 

0.71 ± 0.45 ± 0.04 (< 1 . 00 ) 

0.828 

1.056 

9.0 ± 3.3 ± 0.6 

4.210 

54.6 

21.5 

21.0 

20.8 

20.9 

< 18.9 

26.6 ± 5.4 

0.42 ± 0.21 ± 0.03 (< 0 . 72 ) 

0.813 

1.057 

22.7 ± 4.6 ± 1.5 

4.220 

54.1 

21.6 

21.2 

20.8 

21.1 

< 12.6 

31.9 ± 5.7 

0.18 ± 0.14 ± 0.02 (< 0 . 41 ) 

0.810 

1.057 

27.4 ± 4.9 ± 1.8 

4.230 

1091.7 

22.0 

21.1 

21.0 

21.0 

236.8 ± 25.0 

825.1 ± 29.8 

0.28 ± 0.03 ± 0.02 

0.805 

1.056 

35.4 ± 1.3 ± 2.2 

4.245 

55.6 

22.3 

21.6 

21.1 

21.5 

< 15.2 

49.0 ± 7.1 

0.15 ± 0.10 ± 0.02 (< 0 . 32 ) 

0.806 

1.056 

40.3 ± 5.8 ± 2.7 

4.260 

825.7 

22.6 

21.2 

21.4 

21.2 

73.1 ± 16.5 

507.3 ± 23.4 

0.14 ± 0.03 ± 0.01 

0.815 

1.054 

28.3 ± 1.3 ± 1.8 

4.310 

44.9 

22.5 

20.4 

20.7 

20.5 

< 7.9 

25.5 ± 5.1 

0.07 ± 0.12 ± 0.01 (< 0 . 29 ) 

0.916 

1.052 

24.1 ± 4.9 ± 1.6 

4.360 

539.8 

21.5 

18.8 

19.1 

18.9 

41.8 ± 10.8 

182.8 ± 14.2 

0.20 ± 0.05 ± 0.02 

1.038 

1.051 

13.8 ± 1.1 ± 0.9 

4.390 

55.2 

21.4 

17.7 

18.4 

17.7 

< 5.2 

6.2 ± 2.6 

0.00 ± 1.02 ± 0.00 (< 0 . 71 ) 

1.088 

1.051 

4.7 ± 1.9 ± 0.3 

4.420 

44.7 

21.7 

16.8 

17.9 

16.8 

< 3.8 

2.9 ± 2.1 

0.00 ± 0.56 ± 0.00 (< 1 . 00 ) 

1.132 

1.053 

2.7 ± 1.9 ± 0.2 


The uncertainties in R and a Bom include contributions 
from the luminosity (0% for R and 1.0% for <JBorn) 0; 
tracking efficiency (0% and 2.0%) [l^, 7r° selection ef¬ 
ficiency (0% and 4.0%) [1^, muon identification effi¬ 
ciency (0% and 3.0%), background shape (3.0% and 
0.6%), peaking backgrounds (1.4% and 1.4%), fitting 
range (2.6% and 0.6%), kinematic fit (2.2% and 1.7%), 
intermediate-state branching ratios (0% and 0.5%), sig¬ 
nal parameterization (1.9% and 1.9%), input cross sec¬ 
tion line shape in KKMC (0% and 0.6%) [U, [^, line 
shape of e+e“ —>■ (1.1% to 12.3% and 0% to 3.2%, 

depending on ifcm), and decay models of in the 

MC (0.2% to 6.3% and 0.2% to 6.3%). An uncertainty 
of 0% in R signifies that the effect of that source of sys¬ 
tematic uncertainty cancels in the ratio. Results for R 
and a Born with systematic errors are given in Table ID 
In cases where there is no statistically significant signal, 
upper limits are defined as sums of 90% confidence level 
statistical upper limits plus systematic errors. 


In summary, we have observed a new charmonium- 
like state .Zc(3900)° in e+e” —?► 7r°7r°J/'0 with a sta¬ 
tistical significance of lO.dcr. The mass and width of 
Zc(3900)° are measured to be 3894.8 ± 2.3 ± 3.2 MeV/c^ 
and 29.6 ±8.2 ±8.2 MeV, respectively. We interpret this 
state as the neutral partner of the four-quark state can¬ 
didate Zc(3900)^, since it decays to TT^J/ip and its mass 
is close to the mass of ^c(3900}^. The previous report 
of 3.5(7 evidence for Zc(3900)° 0 included values of the 
mass and width that are consistent with our results, but 
are much less precise. We have also measured the cross 
section ratio R = ^ and the 

(7{e+e 77^ J/yj) 


Born cross section for e“'"e“ —>■ in the energy 

range from 4.190 to 4.420 GeV. The measured Born cross 
sections are about half of those for e+e” —>■ tt+tt” J/^/; 
that were measured by Belle 0 , consistent with the 
isospin symmetry expectation for resonances. 
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